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ABSTRACT 

In this study we present first results from multi-wavelength Hubble Space 
Telescope [HST) observations of the Galactic globular cluster (GG) NGG2808 
as an extension of the Hubble Space Telescope UV Legacy Survey of Galactic GCs 
(GO-13297 and previous proprietary and HST archive data). Our analysis al¬ 
lowed us to disclose a multiple-stellar-population phenomenon inNGG2808 even 
more complex than previously thought. We have separated at least hve different 
populations along the main sequence and the red giant branch (RGB), that we 
name A, B, G, D and E (though an even hner subdivision may be suggested by 
the data). We identihed the RGB bump in four out of the hve RGBs. To explore 
the origin of this complex GMD, we have combined our multi-wavelength HST 
photometry with synthetic spectra, generated by assuming different chemical 
compositions. The comparison of observed colors with synthetic spectra suggests 
that the hve stellar populations have diherent contents of light elements and he¬ 
lium. Specihcally, if we assume that NGG2808 is homogeneous in [Fe/H] (as 
suggested by spectroscopy for Populations B, G, D, E, but lacking for Popula¬ 
tion A) and that population A has a primordial helium abundance, we hnd that 
populations B, G, D, E are enhanced in helium by A Y~0.03, 0.03, 0.08, 0.13, 
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respectively. We obtain similar results by comparing the magnitude of the RGB 
bumps with models. Planned spectroscopic observations will test whether also 
Population A has the same metallicity, or whether its photometric differences 
with Population B can be ascribed to small [Fe/H] and [0/H] differences rather 
than to helium. 

Subject headings: stars: Population II — globular clusters individual: NGC 2808 


1. Introduction 

Recent studies, based on multi-wavelength photometry, have revealed that the color- 
magnitude diagram (GMD) of all Galactic globular clusters (GGs) so far explored (Piotto 
et al. 2015, hereafter Paper I) is made of distinct sequences of stars that can be traced con¬ 
tinuously from the bottom of the main sequence (MS) up to the tip of the red giant branch 
(RGB) and through the horizontal branch (HB) and the asymptotic giant branch (AGB). 
These sequences stand in contrast to the traditional view of GGs as the best example of 
simple stellar populations, i.e. made of stars born all at the same time and with the same 
chemical composition, conhrming previous findings of GN variations in MS and RGB stars 
(Gannon et al. 1998; Grundahl et al. 1998; Grundahl 1999). 

The Hubble Space Telescope UV Legacy Survey of Galactic GCs is an Hubble Space Telescope 
(HST) project to observe 54 GGs through the hlters F275W, F336W, F438W of the Wide 
Field Gamera 3 (WFG3) onboard HST (see Paper I). This dataset complements the existing 
F606W and F814W photometry from the Advanced Gamera for Survey (GO-10775, Sara- 
jedini et al. 2007; Anderson et al. 2008) and is specihcally designed to map multiple stellar 
populations in GGs. 

NGG2808 is one of the most intriguing Galactic GGs in the context of multiple stellar 
populations. It hosts a multimodal MS (D’Antona et al. 2005; Piotto et al. 2007; Milone et 
al. 2012a; and Paper I), and exhibits a multimodal HB (Sosin et al. 1997; Bedin et al. 2000; 
Dalessandro et al. 2011) and RGB (Lee et al. 2009; Monelli et al. 2013; Paper I). Spectroscopy 
has shown star-to-star variations of several light elements, lithium and an extended Na-0 
anticorrelation (Norris & Smith 1983; Garretta et al. 2006, 2010; Gratton et al. 2011; Marino 
et al. 2014; Garretta2014; D’Orazi et al. 2015). These observations have been interpreted 


^Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
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with multiple populations of stars with different helium abundance, from primordial abun¬ 
dance, Y ~0.246, up to extreme enhancement, Y ~0.38 (e.g. D’Antona et ah 2002, 2005; 
Piotto et ah 2007; Milone et ah 2012a). Evidence of helium enhancement in NGC2808 has 
been also confirmed by direct measurements of helium-rich stars along the RGB and the HB 
(Pasquini et ah 2011; Marino et ah 2014). 

While previous studies on multiple MSs in NGG2808 were based on visual or near- 
infrared photometry, in this paper we extend the study to the ultraviolet. The ultraviolet 
region of the spectrum is indeed very powerful in the study of multiple stellar populations 
with different chemical composition. Molecular bands, such as OH, NH, GH, and GN affect 
the ultraviolet and blue wavelengths, that are thus sensitive to populations with different G, 
N, and O compositions (Milone et al. 2012b; Paper I). 

In this paper we use multi-wavelength ultraviolet and visual photometry (from Paper I 
) of stars in a field centered on NGG2808 in order to identify multiple stellar populations in 
the GMDs. The behavior of multiple sequences in appropriate GMDs made with different 
combinations of colors and magnitudes will provide unique information on the helium and 
light-element content of the different stellar populations of this extreme GG. 

The paper is organized as follows. In Sect. [2] we present the data and the data reduction. 
In Sect.[3]we analyze the GMDs and investigate multiple populations along the RGB and the 
MS. Helium and G, N, O abundances of the stellar populations are inferred in Sect. |3]from 
multiple MS and RGB locations in the GMD. The bump of multiple RGBs of NGG2808 are 
analyzed in Sect. [5], while Sects. [7] and [8] are dedicated to the AGB and the HB. A discussion 
will follow in Sect. E 


2. Data and Data Reduction 

In our study of NGG2808 we have used archival and proprietary images taken with the 
Wide Field Ghannel of the Advanced Gamera for Surveys (WFG/AGS) and the Ultraviolet 
and Visual Ghannel of the Wide Field Gamera 3 (UVIS/WFG3) on board of the Hubble 
Space Teleseope {HST). Table 1 summarizes the dataset. 

The poor charge-transfer efficiency (GTE) in the UVIS/WFG3 and AGS/WFG images 
have been corrected by following the recipe by Anderson & Bedin(2010). Photometry and 
astrometry of UVIS/WFG3 images were already presented in Paper I and were obtained 
with img2xymJJVISH9 x 10, which is a software package presented by Bellini et al. (2010) 
and mostly adapted from img2xymJYFI (Anderson et al. 2006). We used pixel-area and 
geometric-distortion corrections from Bellini Sz Bedin (2009) and Bellini, Anderson & Be- 
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din (2011). The photometry has been calibrated as in Bedin et ah (2005), and is using the 
encircled energy and the zero points available at the STScI web page. We used the photo¬ 
metric and astrometric catalogs from WFC/ACS data published by Anderson et ah (2008), 
Sarajedini et ah (2007), Milone et ah (2012a), Piotto et ah (2007), which were obtained from 
GO-9899, GO-10922, and GO-10775 WFG/AGS data. 

In order to investigate multiple stellar populations in NGG2808 we are interested in 
stars for which high-accuracy photometry is available. The stellar catalogs were purged of 
stars that are poorly measured by using the procedure described by Milone et al. (2009) 
and based on the quality indexes provided by our software (see Anderson et al. 2006, 2008). 
Photometry has been corrected for differential reddening following the recipe in Milone et 
al. (2012a). 


INSTR. 

DATE 

NxEXPTIME 

FILTER 

PROGRAM 

PI 

ACS/WFC 

May 5 2004 

6 X 340s 

F475W 

9899 

G. Piotto 

ACS/WFC 

Aug 9 and Nov 2 2006 

20s-|-2x350s-|-2x360s 

F475W 

10922 

G. Piotto 

ACS/WFC 

Aug 9 and Nov 1 2006 

lOs+3 X 350s+3 X 360s 

F814W 

10922 

G. Piotto 

ACS/WFC 

Jan 1 2006 

23sH-5x360s 

F606W 

10775 

A. Sarajedini 

ACS/WFC 

Jan 1 2006 

23sH-5x370s 

F814W 

10775 

A. Sarajedini 

WFC3/UVIS 

Sep 08-09 2013 

12 X 985s 

F275W 

12605 

G. Piotto 

WFC3/UVIS 

Sep 08-09 2013 

6 X 650s 

F336W 

12605 

G. Piotto 

WFC3/UVIS 

Sep 08-09 2013 

6x97s 

F438W 

12605 

G. Piotto 


Table 1: List of the data sets used in this paper. 


3. The multiple photometric components along the CMD of NGC2808 

As already discussed in Sect. [H previous studies based on ACS/HST and ground-based 
photometry have shown that NGG 2808 has at least a triple MS (Piotto et al. 2007; Milone et 
al. 2012a) and a broadened RGB (Lee et al. 2009; Monelli et al. 2013; Paper I). An inspection 
of the large number of GMDs that we derived from six-band photometry immediately reveals 
that NGG 2808 is even more complex than initially thought, and hosts more than three stellar 
populations. 

A visual example of its complexity is provided by the GMDs in Fig. [T] This hgure 
shows several diagrams, after the quality selection and the differential-reddening correction 
described in the previous section were applied. To derive some of the diagrams of Fig. [1] 
we have dehned the pseudo-magnitudes mF 336 w,F 275 W,F 8 i 4 W = ("^Fasew — ^F 275 W + "^fsmw) 
and mF 275 w,F 336 W,F 8 i 4 W = ('"iF 275 W “ ^F 336 W + m-Fsuw): which allow US to better distin¬ 
guish multiple sequences along the RGB and the MS. An inspection of these GMDs im¬ 
mediately suggests that both the RGB and MS are made of multiple sequences, that look 
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discrete in the mF275W vs. mF275W — ’^FSUW, '^F275W,F336W,F814W vs. mF275W — '"2-F336W, and 
mF 336 W,F 275 W,F 8 i 4 W VS. 2 mF 275 w “ nrF 438 W ” n 7 .F 8 i 4 W diagrams. We also observe a widely- 
spread SGB in the mF 275 w vs. mF 336 w ~ 'nrF 438 w CMD as shown in the upper-right inset 
of Fig. [TJ In the following we discuss the observed morphology of the CMD at various 
evolutionary stages: the RGB, MS, AGB, and HB. 


3.1. The quintuple RGB 

Along the RGB, the behavior of multiple populations dramatically changes from one 
CMD to another. In order to investigate this phenomenon, in Fig. [2]we compare the mF 8 i 4 w 
vs. mF 275 w — 'nrF 8 i 4 W CMD of RGB stars (upper-left panel) and of the mF 8 i 4 W vs. mF 336 W — 
nrF 438 W CMD (upper-right panel). The insets show the Hess diagrams for stars in the 
magnitude interval with 14.5 < mF 8 i 4 w < 17-7 where multiple RGBs are clearly visible. 
Multiple stellar populations manifest themselves as four separate sequences in mF 275 w ~ 
mF 8 i 4 W) while the mF 336 w ~ ''^F 438 W color distribution is more broadened and only two 
RGBs can be recognized. 

In order to compare the two CMDs, we used the procedure introduced by Milone et 
al. (2015, hereafter paperII) in their study of multiple populations in M2, and illustrated 
in Fig. [2] for the case of NGC2808. For that purpose, we drew two hducial lines in each 
CMD. The blue and the red fiducial mark the bluest and the reddest envelope of the RGB, 
respectively and have been derived as follows. We have divided the RGB portion with 
^77 f 814W >14.8 into intervals of 0.2 magnitudes in F814W band. For each interval, we have 
determined the 4*^ and the 96*^ percentile of the mF 275 w ~ '^F 814 W or rnF 336 W ~ 
color distribution and the median mF 8 i 4 w magnitude. The points corresponding to the 4*^ 
percentile and the median magnitude have been interpolated with a cubic spline to derive 
the blue fiducial, while the red fiducial has been similarly derived. Due to small number 
statistics it is not possible to infer robust estimates of the RGB envelopes with this method 
at brighter luminosities. Therefore, the portions of the blue and the red line in the magnitude 
interval with mF 8 i 4 w <14.8 have been derived by hand trying to follow the blue and the red 
envelope of the RGB, respectively. 

Then we have verticalized the two CMDs in a way that the blue and the red fiducials 
translate into vertical lines with abscissa —1 and 0, respectively. To do this, we defined for 
each star: 

^x=[(^~^blue fiducial)/(Wred fiducial—Wblue fiducial)]-! where X= (mF275W—'"2-F814w), (mF336W — 

mF 438 w) and Xbiue fiducial and Xred fiducial are obtained by subtracting the color of the hducial 
at the corresponding F814W magnitude from the color of each star. The verticalized mF 8 i 4 w 
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Fig. 1.— 'n^F 275 W vs. mF 275 W ~ '^F 8 i 4 W CMD of NGC2808. The 7nF275w,F336W,F8i4W against 
niF 275 W—'niF 336 W (bottom-left inset), mF336W,F275W,F8i4W against 2 mF 275 W—'niF438W—niFsuw 
(bottom-right inset), and mF 275 W vs. mF 336 W —'niF 438 W (upper-right inset) diagrams highlight 
multiple sequences along the RGB, the MS, and the SGB, respectively. 
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vs. Ap 275 Yv F 814 W mF 8 i 4 W VS. Apggg-y^^ F 438 W diagrams are plotted in the lower-left panels 
of Fig. [2j RGB stars in NGC2808 are clnstered aronnd distinct values of Apggg.^ p 43 g.^ and 
^F 275 W F 814 W’ showu in the bottom-right panel of Fig. [21 

As previously discussed by Anderson et ah (2008, see their Sect. 8.1), F814W photometry 
of bright RGB stars is less accurate than that of the remaining RGB stars because it has 
been derived by using saturated stars (see Anderson et ah 2008 for details). Indeed, multiple 
sequences are less evident above the gray dashed lines in the lower-left panels of Fig. [21 
Dashed lines are placed at mp 8 i 4 w = 14.68. To investigate whether the distinct sequences 
can be also detected along the brightest RGB segment or not, we have marked stars with 
’^F 8 i 4 W < 14.68 with red dots in the lower panels of Fig. [21 The distribution of these bright 
RGB stars on the Apggg.^^ f 438 W ^F 275 W F 814 W shows that these stars share the same 

color distribution as the fainter RGB stars. 

To further investigate the stellar populations along the RGB, in the upper-right panel 
of Fig. [3] we plot the Ap 33 g.yy p 43 g.yy vs. Ap 275 yy pg 44 yy Hess diagram. At least five main clumps 
of RGB stars are clearly visible. These are selected by eye and designated A, B, G, D, 
and E and are colored green, orange, yellow, cyan, and blue, respectively (see the lower-left 
panel diagram). These color codes will be consistently used in the paper. RGB-A-E contain 
5.8±0.5%, 17.4±0.9%, 26.4±1.2%, 31.3±1.3%, and 19.1±1.0% of the total number of RGB 
stars with 12.25 < rripsuw < 17.70, respectively. In Section 13731 we show that populations 
A-E have different chemical composition. 

The Ap 275 yy P 814 W ^^d Ap 33 gyy P 438 W distributions of RGB stars are shown in the upper- 
left and lower-right panel of Fig. [3l respectively. Black histograms represent the whole 
sample of RGB stars shown in the lower-left panel, while the distributions of the five distinct 
RGBs are plotted with shaded-colored histograms. The Ap 275 yy p 8 i 4 'w ^md Apgggyy p 43 gyy dis¬ 
tributions exhibit significant differences. The histogram distribution of Ap 275 w,F 8 i 4 W clearly 
shows three main peaks at Ap 275 yyp 8 i 4 w ~ —0.9,—0.6, and —0.3. The hrst and the second 
clumps are mainly composed of population-E and population-D stars, respectively, while 
the third peak is a mix of both population-B and population-G stars. A less populous 
peak, corresponding to population A, is located at Ap 275 yyp 8 i 4 w ~0.0. In contrast, the 
^F 336 WF 438 W distribution looks bimodal. Most of the stars of populations G, D, and E have 
^F 336 W F 438 W > —0.5 and determine the main peak at Apgggyy p 438 w ~ —0.2. A second peak, 
mostly composed of popuIation-B stars, is located around Apgggyy p 43 gyy ~ —0.8. 

In addition we note that: 

• Populations B and G are mixed in the Ap 275 yy p 8 i 4 'w color range while they have distinct 
'^F 336 W,F 438 W values, with population-B stars having also smaller Apgggyy p 43 gyy values. 
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Fig. 2.— Upper panels: Zoom of the mF 8 i 4 w vs. mF 275 w ~ ''^F 814 W (left) and of the mF 8 i 4 w 
vs. mF 336 w — ''71F438W (right) CMD of NGC2808 around the RGB. Only RGB stars colored 
black are used in the following analysis. Red and blue lines are the hducials adopted to 
verticalize the RGB (see text for details). The insets show the Hess diagram for RGB stars 
with 12.25 < mF 8 i 4 w < 17.7. Lower panels: Verticalized mF 8 i 4 w vs. Ap 275 -\Y f 814 W (left) and 
mF 8 i 4 W vs. Ap 33 g-^ F 438 W (middle) diagrams for RGB stars. Ap 33 g^ f 438 W 1® plotted against 
Ap 275 w F 814 W lower-riglit panel. RGB stars with mF 8 i 4 w < 14.68 are colored red. 
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• Population-A stars have larger Ap 275 Yv f 814 W than both populations B and C. The color 
order is different in Af 336 w,F 438 W; where the histogram of population-A stars is located 
between the histograms of population B and C. 

• Since, the analyzed RGB stars cover the same F814W magnitude interval and have 
similar mF 336 w aiid mF 438 w magnitudes, their photometric errors are similar. We note 
that the A^ 33 g^_F 438 W spread for populations B and C, (cr^F 336 W,F 438 W =0.13±0.01, 
and crAK 36 WF 438 W =0.16±0.01) are significantly larger than the spread observed for 
population-A, D and E stars (crAF 336 W,F 438 W ~0.10i0.02, <^af 336 W,F 438 W ~0.10i0.01, 
and cta^gw F 438 W =0.10±0.01). This fact indicates that the mF 336 W — rnF 438 w color 
spread observed for population B and C is, in part, intrinsic and that both groups B 
and group C are not simple stellar populations. In fact, a visual inspection of the Hess 
diagram of Fig. [3] suggests that both groups consist of two clumps of stars which are 
clustered around Af 336 wf 438 W ~ —0.85, —0.65 (group B) and —0.2, and —0.1 (group 
C) thus suggesting that stars in both groups B and C do not have homogeneous chem¬ 
ical composition. More data are needed to establish whether these clumps correspond 
to distinct stellar populations. 

The causes of the ‘discreteness’ of multiple populations as observed in the CMD and 
two-color diagram of some GCs are still unknown, and have been associated with distinct 
bursts of star formation (see Renzini2008 for a critical discussion). The referee has pointed 
out that the distinct bumps in the diagrams of Fig. [3] could indicate that some abundances are 
favored over the others and suggested a possible connection between the abundances of stars 
in the distinct clumps of NGG2808 and metal mixtures that are consistent with equilibrium 
GN or equilibrium ON cycling. While this hypothesis deserves some investigation that is 
beyond the purposes of our paper, we emphasize that the evidence of discrete populations 
in NGG2808 provides strong constraint for any model of formation and evolution of stellar 
populations in GGs. 

Upon request of the referee and of the Statistical Editor oi the journal. Prof. Eric Feigel- 
son, we have used the Mcluster GRAN package in the public domain R statistical software 
system to estimate how many groups are statistically signihcant. This package is based on 
the method described in details in the monograph ‘Finite Mixture Models’ by McLachlan Sz 
Peel (2000). It performs the maximum likelihood hts to different number of stellar groups, 
and evaluate the number of groups by the Bayesian Information Griterion (BIG) penalized 
likelihood measure for model complexity. 

To do this it uses several different assumptions about shape and size of the different 
populations in a plot such as that shown in Fig. O For each shape and size that we adopted 
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for the populations, we assumed a number, N, of stellar populations from 1 to 20 and 
estimated a BIC for each combination. We obtain the best BIC value (BIC=1784) for N=6 
under the assumption the the stellar populations have equal shapes but variable volume 
and orientations (VEV). The second most-likely explanation (BIC=1778) corresponds to 
N=6 but assumes equal shape, volume, and orientation. The third-best value (BIC=1776) 
corresponds to a VEV assumption and seven stellar populations. All the three best models 
assume ellipsoidal distributions. 

Results from this statistical analysis support the conclusion that our observations of 
NGC2808 are consistent with more than Eve groups of stars, and that group C hosts more 
than one stellar population. The third best BIC value suggests that also the group B is not 
consistent with a simple population. Thus, the statistical analysis confirms what was already 
pretty evident from a pure eye inspection of the plots. In the following, we will study the 
five most-evident stellar populations, A-E. 


3.2. Multiple populations along the MS 

The MS of NGC 2808 exhibits different patterns in CMDs based on different photometric 
bands,in close analogy with what we observe along the RGB. This is shown in the upper 
panels of Fig. 0] where we compare the mF 8 i 4 w vs. mF 275 w ~ '"i-Fsuw (left panel) and the 
^F 8 i 4 W vs. mF 336 w ~ ''^F 438 W (right panel) CMDs of MS stars with 19.6 < mF 8 i 4 w < 20.7. 
The MS looks discrete in mF 275 w — ''^F 814 W with three distinct components, in contrast with 
the mF 336 w — ''^F 438 W color distribution, which looks broadened without any evidence for 
discrete sequences. 

In order to identify the different stellar populations, we have verticalized the MSs 
by following the same recipe as introduced in Sect. 13.11 for the RGB, and using the fidu¬ 
cial lines drawn in the upper-panel CMDs. The mF 8 i 4 w vs. Ap 275 Yv fsuw ?^F 814 W 

vs. Ap 33 g^ F 438 W diagrams are plotted in the lower-left and lower-middle panels, while the 
lower-right panel shows A^ 33 gw p 43 gw against A^ 275 W,F 814 W 

The pseudo color C'f275W,F336W,F438W = (’^F 275 W — '^F 336 w) —(^F 336 W “ ''^F438w) defined 
by Milone et ah (2013) is another valuable tool to identify multiple populations in GCs. To 
better distinguish the distinct MSs and RGBs of NGC 2808 we show in the left panel of 
Fig. O the mF 8 i 4 w vs. Gf 275 w,F 336 W,F 438 W pseudo CMD for this cluster. The red and the blue 
lines superimposed on this diagram are the envelopes of the MS and the RGB, and have 
been determined with the same procedure as in Sect. 13.11 These two fiducials are used to 
verticalize the MS and the RGB. The verticalized mF 8 i 4 w vs. A^p 275 -y\^ p 33 g-\Y p 43 gYv diagram 
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Fig. 3. — Reproduction of the Af336W,f438W vs. Af275w,F814W diagram of Fig. El Stars in the 
A, B, C, D, and E group are colored green, orange, yellow, cyan, and blue, respectively (lower- 
left). The corresponding Hess diagram is plotted in the upper-right panel. The histograms of 
the normalized Af 275 w,F 814 W and Af 336 w,F 438 W distributions for all the analyzed RGB stars 
are plotted in black in the upper-left and lower-right panel, respectively. The shaded colored 
histograms show the distributions for each of the hve populations dehned in the lower-left 
panel. 























- 12 - 


for RGB and MS stars are plotted in the middle panels. An inspection of these hgnres reveals 
that three distinct seqnences are present along the RGB, while only two MSs are visible. In 
the right panels of Fig. E]we plot A^275W,F814W against A^p 275 w,f 336W,F438W ^^^B (npper 

panel) and MS stars (lower panel). 

In order to identify stellar populations along the MS, we exploit the Ap 33 g^p 43 g^ 
vs. Ap 275 Y\f^F 8 i 4 W diagram and the Ap 275 Y\f f 336W,F438 W Af 275 w,F 814W Hess diagram shown 

in the lower-left and upper-right panels of Fig. El respectively. The distribution of MS stars 
in this plane is similar to what observed for the RGB, as better highlighted by the Hess 
diagram of in the upper-right panel Fig. El There are at least four groups of MS stars that 
we denominate B, G, D, and E, and colored orange, yellow, cyan, and blue, in the bottom-left 
and rightmost panels, in close analogy with what was done for the RGB. Golors introduced 
in this hgure will be used consistently hereafter. Noticeably, the separation among the four 
groups is less clear for the MS than in the case of the RGB. This could be due to fact that 
colors of these relatively hot MS stars are less sensitive to light-element variations than the 
RGB. Indeed we have shown in our previous papers that the color difference between multiple 
MSs and RGBs is due, apart from helium, to different strengths of the molecular bands be¬ 
tween the distinct populations of stars (Marino et al. 2008; Milone et al. 2012). In particular, 
the OH band and the GH G-band, which are stronger in the stellar population with the same 
chemical composition as halo held stars of the same metallicity, mainly fall in the F275W 
and the F438W band, respectively, while the NH band, which is weaker in stars of this 
population, mainly affects the F336W magnitude. Population A is not clearly distinguish¬ 
able, even if a stellar overdensity can be recognized in the Hess diagram at (Ap 275 -w f 814 w! 
Af 336W F438 w) ~(0.0;—0.5). We tentatively associate these stars with population A and color 
them green in the lower-left panel diagram. In the middle- and right-upper panels of Fig. El 
we compare the Ap275Yv f 336W F438 W Af 275 w F 814W Hess diagrams for RGB and MS stars, 
while, in the corresponding lower panels, we show the position of Populations A-E in this 
plane. 

In order to further investigate whether MS-A stars either correspond to a distinct stellar 
population or their position in Ap275-\Y f814W ^F336W F438W pl^ne is entirely due to mea¬ 

surement errors, we adopt a procedure introduced by Anderson et al. (2009) and illustrated 
in Fig. [71 In the left panel we show the mp 8 i 4 w against mF 475 w ~ ''^F 814 W GMD from Milone 
et al. (2012a). The red and blue lines superimposed on the GMD are the fiducials of the red 
and the blue MSs and are drawn by hand. The verticalized mpsuw vs. Ap 475 ^ f 814 W diagram 
is plotted in the central panel, while the right panel shows Ap 275 YVP 8 i 4 w vs. Ap 475 .y^ pg 44 .y^. 
Stars in common with this paper are marked with colored circles. 

Photometry by Milone et al. (2012a) comes from AGS/WFG images, hence represents a 
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Fig. 4.— Upper panels: mpsuw vs. mF 275 w — '"J-fsmw (left) and of the mpsMW vs. mpasew — 
mF 438 w (right) CMD for MS stars in NGC2808. The fiducials nsed to verticalize the RGB 
are represented with red and blue lines (see text for details). Lower panels: Verticalized 
"iF8i4W vs. A^275W,F814W ( 1 ^^) ^ud mpsMW VS. p 43 g^ (middle) diagram for the stars 

in the upper panels. Apggg-^^^ f 438 W 1® plotted against Ap 275 Yv f 814 W 1 ^ lower-right panel. 
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Fig. 5.— Left: mpsuw vs. Cf 275 w,F 336 W,F 438 W diagram for NGC2808. The red and blue 
lines superimposed on the diagram are the hducial lines used to verticalize the MS and the 
RGB. See text for details. Middle: Verticalized MS (upper panel) and RGB (lower panel). 
Right: f 336W F438W ^F275W F814W RGB and MS stars plotted in the middle 

panels. 
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Fig. 6 . — Left: Reproduction of the Af336w,F438W vs. Af275w,F814W diagram of Fig. [ 5 l Stars 
in the A, B, C, D, and E groups, dehned in this hgure, are colored green, orange, yellow, cyan, 
and blue, respectively (lower panel). The corresponding Hess diagram is shown in the upper- 
left panel. Central and middle lower panels show Ap275-\Y f336W F438W against Ap275-\Y f814W 
for RGB and MS stars, respectively while the corresponding Hess diagrams are plotted in 
the upper panels. 
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different dataset than the WFC3 ones nsed in this paper. If the large Ap 275 Yv f 814 W valne for 
stars in the group A derived from UVIS/WFC3 photometry is entirely due to photometric 
errors, then these stars have the same probability of having either small or large Ap 475 ^ f 814 W 
derived from WFC/ACS. The systematically large Ap 475 -^ pg 44 -y\^ value of stars in the group 
A, shown in the right panel of Fig. [TJ is evidence that they belong to a distinct stellar 
population. 



Fig. 7.— Left panel: ruFsuw against mF 475 w — ''^F 814 W from the ACS/WFC photometry 
published by Milone et ah (2012a). Red and blue lines are the hducials of the red and blue 
MS. Middle panel: verticalized mp 8 i 4 w vs. Ap 475 ^ pg 44 ^ diagram. Stars in common with 
this paper are represented with colored circles. Right panel: F814W (from this paper) 

vs. Ap 475 ^ pg 44 ^ (from Milone et al. 2012a). 


3.3. Chemical composition of stellar populations 

Spectroscopy of RGB stars has revealed that NGC 2808 exhibits a very extended sodium- 
oxygen anticorrelation, with [0/Fe] spanning a range of more than 1 dex (Garretta et 
al. 2006). Twenty-seven stars analyzed by Garretta and collaborators are also included in our 
photometric sample thus providing useful information on the chemical composition of the 
stellar populations we have identified in the previous sections. Garretta’s stars are marked 
with large symbols in Fig. |8] where we reproduce the Apggg-^^^ f 438 W vs. Ap 275 -\Y f 814 W diagram 
of Fig. [3] (left panel), the Na-0 anticorrelation from Garretta et al. (2006, middle panel), 
and the Mg-Al anticorrelation from Garretta (2014, right panel). Noticeably, stars in the B, 
G, D, and E stellar groups dehned in this paper have almost the same iron content within 
~0.05 dex but populate different regions of the Na-0 plane. The average elemental abun¬ 
dance for stars of population B, G, D, and E is listed in Table [21 where we also provide the 
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dispersion, a, and the number of stars, N, in each population with available abundances. 
The error is estimated as the a divided by the square root of N—1. Unfortunately, there are 
no population-A stars in the sample analyzed by Carretta and collaborators. 

The hve population-B stars in common with the Carretta et ah sample have all pri¬ 
mordial Na and O ([Na/Fe]~0.0, [O/Fe]~0.3), and the ten population-C stars are slightly 
enhanced in sodium ([Na/Fe]~0.2 dex) and depleted in O ([O/Fe]~0.15 dex) with respect 
to population B. Population-E stars all have very low oxygen abundance and high [Na/Fe], 
while population-D stars are consistent with an intermediate chemical composition. 

More recently, Carretta (2014) has determined Mg and A1 abundances for 31 RGB stars 
in NGC2808 from UVES spectra and detected a very extended Mg-Al anticorrelation with 
three distinct groups of stars with different content of magnesium and aluminum. There are 
5 stars in common with Carretta (2014), all Mg-rich and Al-poor and with similar [Mg/Fe] 
ratio. On average, the two population-C stars are slightly more Al-rich than the three 
population-B stars by ~0.2 dex, but a larger sample is needed to establish the significance 
of this difference. 

Carretta et ah (2010) have dehned three groups of stars in NGC 2808 on the basis of their 
Na and O abundance. A ‘primordial’ component containing all stars with [Na/Fe]<~0.19, an 
‘extreme’ one with [O/Na]<0.9, and an ‘intermediate’ component with intermediate values of 
Na and O. Our Endings show that the component that they designated ‘primordial’ contains 
at least two difierent stellar populations, i.e. B and G, thus suggesting that NGG2808 has 
experienced a very complex star-formation history. 


4. Carbon, Nitrogen, Oxygen, and Helinm of stellar popnlations from 

multiple MSs and RGBs. 

In order to infer the helium abundance of each stellar population of NGG2808 we 
adapted the method introduced by Milone et al. (2012b) in their study of 47Tucanae to 
the case of NGG2808. This method is based on the comparison between the observed 
colors of the multiple sequences and the predictions from appropriate isochrones and the 
relative synthetic spectra. In this section we will exploit the same technique to estimate the 
abundance of G, N, and O for the five populations of NGG2808. 

In Sects. ITT] and 13^ we have used two-color diagrams from appropriate combinations of 
F275W, F336W, F438W, and F814W magnitudes to identify five stellar populations along 
the RGB and the MS. For simplicity, in the following we will indicate as MS-A—E and 
RGB-A—E the groups of MS and RGB stars of populations A—E. To analyze the behavior 
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[0/Fe] 


[Mg/Fe] 


Fig. 8.— Left panel: reproduction of the p 438 w vs. A^ 275 W,F 814 W diagram of Fig. El 

Stars for which spectroscopic measurement are available are marked with large symbols. 
Middle panel: sodium-oxygen anticorrelation for RGB stars of NGC2808 from Garretta et 
al. (2006). Large orange, yellow, cyan, and blue dots indicate spectroscopic targets of popu¬ 
lation B, G, D, and E, respectively. Stars for which only [Na/Fe] measurements are available 
have been arbitrarily plotted at [0/Fe]=l.l. No population-A stars are present in the Gar¬ 
retta et al. sample. Right panel: magnesium-aluminum anticorrelation from Garretta (2014). 
Population-B and population-G stars are indicated with orange and yellow triangles, respec¬ 
tively. 






Pop. 

abundance 

cr 

# 

[O/Fe] 

A 

— 

— 

0 

B 

0.30±0.05 

0.11 

5 

C 

0.16±0.04 

0.11 

9 

D 

-0.37±0.02 

0.02 

3 

E 

-0.66±0.09 

0.12 

3 

[Mg/Fe] 

A 

— 

— 

0 

B 

0.38±0.02 

0.03 

3 

C 

0.35±0.12 

0.12 

2 

D 

— 

— 

0 

E 

— 

— 

0 

[Al/Fe] 

A 

— 

— 

0 

B 

0.05±0.02 

0.03 

3 

C 

0.23±0.13 

0.13 

2 

D 

— 

— 

0 

E 

— 

— 

0 

[Na/Fe] 

A 

— 

— 

0 

B 

0.03±0.06 

0.13 

5 

C 

0.21±0.03 

0.09 

10 

D 

0.40±0.04 

0.09 

6 

E 

0.79±0.08 

0.18 

6 

[Fe/H] 

A 

— 

— 

0 

B 

-1.13±0.02 

0.04 

5 

C 

-1.08±0.02 

0.06 

10 

D 

-1.12±0.02 

0.04 

6 

E 

-1.10±0.03 

0.08 

6 


Table 2: Average abundance of stars in the five stellar populations of NGC2808 defined in 
this paper. Results are inferred by matching photometry with high-resolution-spectroscopy 
measurements by Carretta et ah (2006) and Carretta (2014). 
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of multiple sequences, we have plotted mF 8 i 4 W against mx — mpsuw-, where X=F275W, 
F336W, F438W, FdTSwE, and F606W, and determined the MS and the RGB hducial lines 
of each population. 

Results are shown in Fig. [9] where we have plotted the RGB and the MS hducials in 
the upper and lower panels, respectively. MS-E exhibits the bluest mx — mpsuw color in 
all the GMDs of Fig. [9l Population D has the second bluest MS, while MS-A is redder than 
the other MSs in all plots. MS-B and MS-G are placed between MS-A and MS-D and their 
relative position changes from one GMD to another. MS-G is slightly bluer than MS-B in 
’^F 275 W —’^F 8 i 4 W and mF 438 W —'^F 8 i 4 W colors but is redder than MS-B in mpasew —''^F 814 W- 
MS-B is almost superimposed on the MS-G in the mF 475 w — ''^F 814 W and mpeoew — ’^fsuw 
colors. Lower panels of Fig. [9] show that the mF 275 w “ '" 2 -F 814 W; ''^F 438 W ~ ’^fsuW; and 
’^Feoew “ '^F 8 i 4 W color order of multiple RGBs and MSs are similar. In the mF 8 i 4 w against 
^F 336 W —'^F 8 i 4 W GMD, RGB-E and RGB-B share the bluest colors, while RGB-A, RGB-G, 
and RGB-D dehne a red sequence. 

In order to further investigate multiple MSs and RGBs, we have calculated the mx — 
"iF 8 i 4 W color difference between each MS (or RGB) hducial and MS-B (or RGB-B) hducial at 
a reference magnitude that we indicate as A(mx—mF 8 i 4 w)- Figure fTOl shows A(mx — 

’^F 8 i 4 w) as a function of the central wavelength of the X hlter for MS (left panel, mp^j^^^v = 
20.25) and RGB hducials (right panel, mp^j^^^ = 16.25). We repeated this procedure for 
mp^]^^w=19-80) 19.95, 20.10, 20.40, and 20.55 for the MS and for mp^^^^^ =15.25, 15.75, 
16.75, and 17.25 for the RGB. We hnd that the color separation between populations A 
and B increases with the color baseline. In the case of both populations D and E, A(mx — 
^F 8 i 4 w) grows monotonically for X=F606W, X=F475W, and X=F438W, then it drops for 
X=F336W and reaches its maximum for X=F275W. Population G exhibits almost the same 
color as population B apart from the case of X=F336W, where population G has negative 
A(mx — ?7lF814w)- 

It has been shown that variations in the abundance of both helium and some light 
elements (e.g. G, N, O) are responsible for multiple MSs and RGBs in GGs (e.g. Marino et 
ah 2008; Sbordone et ah 2011; Milone et ah 2012b; Dotter et al. 2015). In a few GGs multiple 
sequences are also due to iron variations (e.g. Pancino et al. 2000; Piotto et al. 2005; Marino 
et al. 2012; Paper II). To investigate the effect of G, N, O, iron and helium on the GMD of 
NGG2808 we started to compute synthetic spectra that we used as reference for population- 
B stars (from here on ‘reference spectrum’). We used BaSTI isochrones (Pietrinferni et 


^Due to the small number of RGB stars for which F475W is available we have not used this filter for the 
study of multiple RGBs 
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Fig. 9.— The green, orange, yellow, cyan, and blue lines are the MS (upper panels) and RGB 
(lower panels) fiducial lines for the population A, B, C, D, and E, respectively in the mpsuw 
vs. mx — nrF 8 i 4 W plane (X=F275W, F336W, F438W, F475W, and F606W). Horizontal gray 
lines mark the magnitudes at which we have calculated the color distance among the MSs 
and the RGBs. We have not used the F475W filter for the RGB due to the small number of 
stars for which we could measure F475W magnitudes. 
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Fig. 10.— Left panel: A(mx — J^fsuw) color distance between MS-B and MS-A, MS-C, 
MS-D, MS-E (yellow, green, cyan, and bine dots) at =20.25 as a fnnction of the 

central wavelength of filter X. Right panel: Color distance between RGB-B and RGB-A, 
RGB-C, RGB-D, and RGB-E measnred at =16.25 vs. the central A of the X filter. 

ah 2004, 2009) to estimate the effective temperatnre (Teff=5619 K) and the snrface gravity 
(log 5^=4.56) at mpgMw = ^ with helinm (Y = 0.278). We assumed for 

NGC2808 a metallicity of [Fe/H]=—1.14 (Harris 1996, 2010 edition), an average reddening 
E(i? — 1/)=0.19 (Bedin et al. 2000) and a distance modulus (m — M)v=15.59 (Harris 1996, 
2010 edition). We computed a synthetic spectrum with [O/Fe]=0.3, as derived for population 
B in Sect. 13.31 from Carretta et al. (2006) measurements, and [C/Fe]=—0.3, and [N/Fe]=0.5 as 
inferred by Bragaglia et al. (2010) for one red-MS star analyzed in their paper. We assumed 
Y=0.278 for population E§. 

We have also computed additional synthetic spectra for stars with the same F814W 
magnitude, but different chemistry. We denominated these spectra as ‘comparison spectra’. 
We assumed for each of them the same chemical composition as the reference spectrum, but 
different abundances for He, C, N, O, and Fe. 


•^This choice for the helium content of population B is somehow arbitrary, and in practice adopted to avoid 
a helium content formally smaller than the primordial one for population A. We note that, in this section, we 
measure helium differences among the different populations, not absolute helium values. A different choice 
of primordial helium content for the reference population would have a negligible impact on the estimated 
AF. 
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Since stars with the same luminosity in the F814W band, but different He or [Fe/H], have 
also different temperature and gravity (e.g., Sbordone et ah 2011; Cassisi et ah 2013), when 
the content of helium and iron are varied, atmospheric parameters are changed accordingly, 
as predicted by BaSTI isochrones. Specihcally, in order to determine the appropriate value 
of Teff = Tgg for a star with magnitude m-psuw = ''^f 814 W ^ given content of helium 

Y = Y*, we have calculated the values of Teg isochrone at mF 8 i 4 W = '^f 814 W from each available 
isochrone and estimated by linearly interpolating the value of Teff^isochrone that corresponds 
to H = Y*. Interpolations are done within the various available isochrones with helium 
frisochrone=0.246, 0.28, 0.30, 0.33, 0.35, 0.40. 

The adopted range of light elements and iron matches observations from high-resolution 
spectroscopy (Carretta et ah 2006; Bragaglia et ah 2010; Gratton et al. 2011; Marino et 
al. 2014). The difference in temperature, and chemical composition between the comparison 
and the reference spectrum are indicated in the left panels of Fig. [TT] 

We used ATLAS12 (Kurucz2005, Castelli 2005, Sbordone et al. 2007) and SYNTHE 
codes (Kurucz2005) to account for the adopted chemical composition and to perform the 
spectral synthesis in the wavelength interval between 2,000 and 10,000 A. Synthetic spectra 
have been integrated over the transmission curves of the F275W, F336W, F438W hlters of 
UVIS/WFC3 and the F475W, F606W, F814W hlters of ACS/WFC, and, for each spectrum, 
we calculated the color mx — wrF 8 i 4 W- 

Results are illustrated in the right panels of Fig. [TT] and provide an indication on the ef¬ 
fect of varying C, N, O, Fe, and He on the observed colors. Left panels show the hux ratio be¬ 
tween the comparison spectrum and the reference spectrum of MS stars with =20.25. 

The corresponding A(mx — mF 8 i 4 w) color differences are plotted in the right panels. 

In the top-left left panels of Fig.[TT]we analyze two comparison spectra sharing the same 
atmospheric parameters and chemical composition as the reference spectrum, but depleted in 
carbon. A difference in [C/Fe] of 0.4, as observed for red- and blue-MS stars by Bragaglia et 
al. (2010), mainly affects mF 438 W — ’^fsuw and marginally changes the other colors studied 
in this work. On the other hand, variations in [N/Fe] produce a large mpssgw ~ '" 2 -F 814 W 
difference (panel b). Because of the effect of the OH, NH, and CH bands on the different 
photometric bands discussed in Sect. 13.21 changing [0/Fe] mainly affects mF 275 W — '^F 814 W 
and mF 336 W — ''^F 814 W colors, with negligible changes for visual colors (panels c). 

In panels d we changed C, N, and O in a way that the overall C-l-N-l-0 abundance in 
the reference and in the comparison spectra remains unchanged. In this case we have large 
differences in 'm.F 275 W ~ '^F 814 W and 'U7 .f 336W — 'n7.F8i4W- 

As shown by Sbordone et al. (2011), helium variation marginally affects the atmospheric 
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structure and the resulting flux distribution. However, stars with the same luminosity but 
different helium have different effective temperature (see Fig. 1 by Sbordone et ah 2011, 
as an example). For this reason, helium variation results in a difference of both visual and 
ultraviolet colors as shown in panels e, with the mF 275 w ~''^F 814 W aiid mpssew ~’^F 814 W color 
differences being significantly larger than those in mF 438 W — ''^F 814 W and mF 475 w — ''^F 814 W- 
Also variations in [Fe/H] affect all the colors studied in this paper, but, in this case, the 
differences corresponding to visual and ultraviolet colors are similar (panels f). 

We conclude that the F275W and F336W bandpasses are very sensitive to the detailed 
chemical abundance of the cluster stars, therefore maximizing the separation among the 
various sub-sequences due to the variations of light-elements and helium. Specifically, the 
F275W band is mainly affected by oxygen variations via the strength of the OH molecular 
bands, while F336W is mostly sensitive to nitrogen via the the NH molecular bands. In 
contrast, the optical colors ('m,F 438 w — ^^F 814 W) ^^F 475 W — ^^F 814 W) and ruFeoew — ^^F 814 w) 
are less affected by C, N, O variations, but are very sensitive to temperature variations 
associated with differences in helium. In addition, a difference in [Fe/H] of less than 0.1 
dex, as inferred for NGC2808 from high-resolution spectroscopy by Carretta et ah (2006), 
corresponds to small color variations. 

Having demonstrated that optical and ultraviolet colors are very sensitive to helium 
and light-element variations, we have estimated the chemical composition of population E. 
To do this we have adopted the procedure described above to calculate synthetic spectra 
with different helium abundances, with Y ranging from 0.246 to 0.400 in steps of AY=0.001, 
[C/Fe], [N/Fe], and [0/Fe] from —2.0 to 2.0 dex in steps of A[C/Fe]=0.1; A[N/Fe]=0.1, and 
A[O/Fe]=0.1. The values of Tgs and logg corresponding to different Y have been derived 
from isochrones as described above. We used chi-square minimization to determine the best 
fit between the synthetic colors and observations. The helium, carbon, nitrogen, oxygen 
difference corresponding to the best-fit model are listed in Table [3] for different values of 
together with its effective temperature and gravity. 

Multi-wavelength photometry of RGB stars also is very sensitive to G, N, O and helium. 
To better constrain the chemical composition of the five stellar populations of NGG2808, 
we have extended the method above described for MS stars to the RGB. We note that the 
values of [G/Fe] and [N/Fe] adopted for population-B stars come from spectroscopy of MS 
stars. When a star ascend the RGB it is affected by mixing phenomena which alter the 
original surface abundance of G and N, while keeping constant G-l-N (e.g. Iben 1967). To 
account for this phenomenon, we followed the recipe by Milone (2015) and assumed that the 
reference spectrum for a RGB-B star has a 0.3-dex lower G abundance than that inferred from 
the MS by Bragaglia et ah (2010). Nitrogen abundance has been determined by assuming 
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Fig. 11.— Left panels: the ratio between the flux of some comparison spectra and the refer¬ 
ence spectrum for a MS-B star at =20.25 is plotted as a function of the wavelength. 

Each comparison spectrum has the same chemical composition as the reference spectrum 
apart from the abundance of one element, as indicated in each panel. When Y and [Fe/H] 
are changed also the atmospheric parameters of the comparison spectrum are varied accord¬ 
ingly, as indicated. Right panels: color difference (A(mx —'mF 8 i 4 w)) between the comparison 
spectra and the reference spectrum (see text for details). 
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that the sum C+N remains constant and corresponds to [N/Fe]=0.62 dex. The adopted 
variation approximately matches the value predicted by Angelou et ah (2011) for multiple 
stellar populations in M 3. 

Results for populations A, C, D, and E are illustrated in Fig. [121 Upper-left panel 
shows the flux ratio between the best-fitting MS-E-comparison spectrum and the reference 
spectrum for =20.25, while the normalized throughput of the filters used in this paper 

is plotted in the middle-left panel. In the lower-left panel we overplotted on the observed 
color difference between MS-B and MS-E of Fig. [10] (blue dots) the color differences derived 
from synthetic spectra. Right panels compare the observed MS and RGB color distance 
between population B and populations A, C, D, and E of Fig. [10] and the color differences 
from the best-£t models of Table [3] and [H 

From the analysis of multiple MSs, we can infer that populations D and E are highly 
helium enhanced by AY~0.11 and AY~0.06, respectively, and are both strongly enhanced 
in N and depleted in O. Population C shares almost the same helium as population B 
(AY~0.01) and is more nitrogen rich and more oxygen poor than the latter. Under the 
assumption of populations A and B having the same metallicity, population A would have 
a lower helium content than population B, by AY ~ 0.03, while it would have almost the 
same nitrogen. 

As suggested by the referee, we have also repeated the helium abundance estimates 
using only the mpeoew ~ '^F 814 W color difference between each sequence and isochrones 
with different helium abundance, as this color is not significantly affected by light-element 
variations (Sbordone et ah 2011). The results are fully consistent with what found from the 
analysis above, within helium differences less than 0.005. 

Results on the estimated compositions of the different populations are shown in Fig. [T31 
where we used colored circles to illustrate the differences in carbon (left panel), nitrogen 
(middle panel), and oxygen (right panel) as a function of the helium variation between 
populations A, C, D, E and population B as inferred from multi-wavelength photometry of 
multiple MSs (upper panels) and multiple RGBs (lower panels). Overall, we see a correlation 
between nitrogen and helium, while both [G/Fe] and [0/Fe] anticorrelate with Y. Noticeably 
there is a large difference of nitrogen between population B and G associated with a small 
helium variation. 

In the right panel of Fig. [T3] colored triangles are obtained from the average oxygen 
abundance determined in Sect. 13.31 from high-resolution spectroscopy of RGB stars by Gar- 
retta et al. (2006) for populations B, G, D, and E. In general, results from spectroscopy and 
photometry are in agreement at one-sigma level, although we note that photometry from 
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Fig. 12.— Left Panels: Flux ratio between the best-fit comparison spectra of population E 
and population B for MS stars at =20.25 is plotted as a function of the wavelength 

in the upper panel where we also show the normalized transmission curves of the hlters used 
in this paper. Observed A(mx — 'mF 8 i 4 w) against the central wavelength of the hlter X are 
plotted with blue circle in the bottom-left panel. Asterisks are inferred from the best-£t 
synthetic spectra. Right Panels Reproduction of Fig. [TOl where we plotted A(mx — ?^f 814 w) 
as a function of the central wavelength of hlter X as observed at =20.25 (top) and 

^m 4 W =16.25 (bottom). In addition, we show the synthetic colors corresponding to the 
best-£t model (see text for details). 






















this paper systematically predicts slightly smaller [0/Fe] variations for populations C, D, 
and E. 
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0.0 

0.0 

C 

5526 

4.59 

0.005 

-0.6 

0.5 

-0.1 

D 

5623 

4.58 

0.051 

-0.8 

0.9 

-0.7 

E 

5747 

4.58 

0.110 

-1.3 

1.2 

-1.1 

^^''F814W“ 

20.55 






A 

5335 

4.61 

-0.033 

0.0 

-0.1 

0.2 

B 

5406 

4.61 

0.000 

0.0 

0.0 

0.0 

C 

5412 

4.61 

0.003 

-0.7 

0.6 

-0.1 

D 

5526 

4.61 

0.056 

-0.8 

0.8 

-0.5 

E 

5655 

4.61 

0.116 

-1.3 

1.1 

-1.0 

AVERAGE 







Pop. 


AY 

A[C/Fe] 

A[N/Fe] 

A[0/Fe] 

A[(C+N+0)/Fe] 

A 


-0.032±0.003 

O.OibO.l 

-O.libO.l 

O.libO.l 

O.libO.2 

B 


0.000 

0.0 

0.0 

0.0 

0.0 

C 


0.002±0.002 

-0.6±0.1 

0.5±0.1 

-O.libO.l 

O.libO.l 

D 


0.051±0.005 

-0.6±0.2 

0.7±0.1 

-0.5ib0.2 

O.OibO.l 

E 


0.106±0.008 

-1.1±0.3 

1.1±0.2 

-0.9ib0.3 

0.3ib0.2 


Table 3: Best-fit atmospheric parameters and relative abundances for synthetic spectra of 
MS stars at different values of We assumed for the MS-B; Y=0.278, [C/Fe]=—0.3, 

[N/Fe]=0.5, [O/Fe]=0.3. See text for details. 
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- 0.05 0 0.05 0.1 - 0.05 0 0.05 0.1 - 0.05 0 0.05 0.1 

AY AY AY 

Fig. 13.— Variation of carbon (left), nitrogen (middle), and oxygen (right) as a function 
of the helium variation for the five populations of NGC2808. Elemental variations are 
calculated with respect to average abundance of population-B stars. Results obtained from 
multiple MSs and RGBs are plotted in the upper and lower panels, respectively. Filled circles 
indicate the He, G, N, O relative abundances inferred from multi-wavelength photometry, 
while the values of A[0/Fe] inferred from spectroscopy of RGB stars are represented with 
triangles. For clarity, spectroscopic measurements are shifted by AY=0.01. 


5. RGB bump of the stellar populations 

In a simple stellar population the luminosity of the RGB bump is an indicator of its 
metallicity, age, helium abundance, and G-l-N-l-0 content (e.g. Gassisi & Salarisl997; Bono 
et ah 2001; Bjork & Ghaboyer 2006; Pietrinferni et ah 2009, Di Gecco et al. 2010; Gassisi 
et al. 2011). In a GMD with multiple sequences as in a GG, the distinct RGB bumps can 
provide information on the age and the chemistry of the various subpopulations. In the 
specihc case of NGG2808, Nataf et al. (2013) found that this GG hosts a broadened RGB 
bump. Although their photometry was not corrected for differential reddening, and therefore 
they did not distinguish the different bumps, they suggested that this peculiar shape of the 
RGB bump is likely due to multiple stellar populations with different helium content. In 
this section we investigate the bump of the RGBs identihed in Sect. 13.11 and infer helium 
abundance by comparing their observed luminosity with theoretical predictions. 
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Pop. 

Tff 

logff 

AY 

A[C/Fe] 

A[N/Fe] 

A[0/Fe] 

"''F814W“ 

15.25 






A 

4905 

2.46 

-0.031 

0.4 

0.1 

0.1 

B 

4933 

2.44 

0.000 

0.0 

0.0 

0.0 

C 

4939 

2.44 

0.007 

0.0 

0.7 

0.1 

D 

4963 

2.43 

0.033 

-0.9 

1.1 

-0.6 

E 

5000 

2.41 

0.074 

-1.0 

1.4 

-0.8 

"''F814W“ 

15.75 






A 

5008 

2.70 

-0.035 

0.0 

0.0 

0.2 

B 

5040 

2.68 

0.000 

0.0 

0.0 

0.0 

C 

5045 

2.68 

0.005 

-0.1 

0.8 

0.1 

D 

5076 

2.66 

0.039 

-0.8 

1.0 

-0.4 

E 

5120 

2.63 

0.086 

-1.1 

1.2 

-0.7 

.v.CUT Z 

^‘•F814W“ 

16.25 






A 

5094 

2.93 

-0.038 

0.3 

-0.1 

0.3 

B 

5132 

2.91 

0.000 

0.0 

0.0 

0.0 

C 

5130 

2.91 

-0.002 

-0.2 

0.7 

0.0 

D 

5175 

2.89 

0.043 

-1.0 

1.1 

-0.4 

E 

5232 

2.86 

0.100 

-1.0 

1.3 

-0.6 

^‘•F814W 

16.75 






A 

5167 

3.14 

-0.034 

0.1 

-0.1 

0.3 

B 

5205 

3.13 

0.000 

0.0 

0.0 

0.0 

C 

5207 

3.13 

0.002 

-0.2 

0.4 

0.0 

D 

5250 

3.11 

0.041 

-0.7 

1.0 

-0.5 

E 

5307 

3.09 

0.092 

-0.7 

1.1 

-0.7 

I-viGUT' — 

"''F814W 

17.25 






A 

5219 

3.36 

-0.037 

0.0 

-0.1 

0.2 

B 

5265 

3.35 

0.000 

0.0 

0.0 

0.0 

C 

5264 

3.35 

-0.001 

-0.2 

0.6 

-0.1 

D 

5316 

3.33 

0.041 

-0.4 

1.1 

-0.7 

E 

5377 

3.31 

0.090 

-0.5 

1.0 

-0.8 

AVERAGE 







Pop. 


AY 

A[C/Fe] 

A[N/Fe] 

A[0/Fe] 

A[(C+N+0)/Fe] 

A 

— 

-0.035±0.003 

0.2±0.2 

O.OitO.l 

0.2±0.1 

0.0±0.2 

B 

— 

0.000 

0.0 

0.0 

0.0 

0.0 

C 

— 

0.002±0.004 

-O.litO.l 

0.6±0.2 

O.OitO.l 

0.2±0.2 

D 

— 

0.040±0.005 

-0.8±0.3 

l.litO.l 

-0.5±0.1 

0.5±0.1 

E 

— 

0.089±0.010 

-0.9±0.3 

1.2±0.2 

-0.7±0.1 

0.5±0.2 


Table 4: Best-fit atmospheric parameters and relative abundances for synthetic spectra 
of RGB stars at different values of We assumed for the RGB-B; Y=0.278, 

[G/Fe] =—0.6, [N/Fe]=0.62, [O/Fe]=0.3. See text for details. 
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The mF 8 i 4 W against mF 275 W—■'^F 814 W and the mF 8 i 4 W against mFsaew—'^F 438 W Hess dia¬ 
grams of RGB stars plotted in Fig. [TTIshow that the RGB bnmp of NGG 2808 exhibits a com¬ 
plex strnctnre. The Hess diagrams reveal mnltiple bnmps, with different F814W Inminosities 
that are associated with the different stellar popnlations. Lower-left and lower-middle panels 
of Fig.mireprodnce the same mFSMW vs. A^ 275 W,F 814 W mpsuw vs. p 438 w diagrams 

of Fig. 131 bnt zoomed-in aronnd the RGB bnmp. In the lower-right panel we plot the his¬ 
togram of the distribntion in mF 8 i 4 w for all the RGB stars shown in the left and in the middle 
panels (black histogram). Golored histograms represent the mF 8 i 4 w distribntion of stars for 
all the stellar popnlations of NGG 2808 we have identihed with the exception of popnlation 
A, which inclndes too few members for an accnrate identihcation of the bnmp magnitnde. 
Histograms have been obtained by adapting the naive estimator (Silverman 1986). Briefly, 
we have hrst dehned a regnlar sample of points (mpg 44 ^), ranging from 14.7 to 15.4, with 
”^m 4 W ~ rnpg 44 ^ = 0.2 mag. Then we have extracted the histogram by associating with 
each point, the nnmber of RGB stars with — u/2 < mFsuw < rnF 8 i 4 w+'^/ 2 - 

We assnmed u =0.013 mag. 

In order to obtain the mp 8 i 4 w Inminosity of the bnmps, we nsed the following procednre. 
We started to estimate by eye a raw position for the bnmp Inminosity in each histogram 
and selected all the points within 0.15 F814W magnitndes from this position. We least 
sqnare htted these points with a Ganssian and considered the position of the center of the 
best-htting Ganssian (mF 8 i 4 w,bump) as a new estimate for the bnmp Inminosity. Finally, we 
selected all the points of the histogram within 2cr (where a is the width of the best-htting 
Ganssian) from mF 8 i 4 w,bump, for a new least-sqnares Ganssian £t. This later determination of 
mF 8 i 4 W,bump corresponds to onr best estimate of the mF 8 i 4 W Inminosity of the bnmp. Errors 
on mF 8 i 4 w,bump are determined by bootstrapping with replacements performed 1,000 times 
on the sample of analyzed RGB stars. The error bars indicate the la (68.27*^) percentile of 
the bootstrapped measnrements. Golored dots in the left and middle panel of Fig. [TTImark 
the position of the bnmps for the four stellar populations in the mF 8 i 4 w vs. Ap 275 w F 814 W 
mF 8 i 4 W vs. Ap 33 g^ F 438 W diagrams. We repeated the same procedures to estimate the bump 
position in the F275W, F336W, F438W, and F814W bands. The RGB-bump magnitudes 
(m^ bump) X=F275W, F336W, F438W, F606W, and F814W) for four stellar populations are 
listed in Tabled while in Fig. [13] we show the difference between the RGB-bump magnitude 
of populations G, D, and E and the RGB-bump magnitude of population B against the 
central wavelength of the X hlter. 

In four bands, namely F275W, F438W, F606W, and F814W, the RGB-bump brightness 
anti-correlates with Ap 275 -\Y p 844 -\Y (e.g. lower left panel of Fig. [T3]). Population E hosts the 
brightest bump. The RGB bumps of populations B and G share almost the same luminosity, 
population-D bump is brighter than those of populations B and G. The magnitude separation 
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Fig. 14.— Upper panels: mpsuw vs. mF275W —''^fsuw (left) and mF8i4W vs. mpsaew—'"iF438W 
Hess diagrams for stars around the RGB bump. Red iso-density contours are superimposed 
on each diagram. Lower panels: zoom of the mF 8 i 4 W vs. Ap 275 w F 814 W (left) and mF 8 i 4 w 
vs. Ap 33 g^ F 438 W (middle) diagram around the RGB bump. The histograms of the F814W 
magnitude distribution are plotted in the right panel for all the stars (black histogram) and 
for the four RGBs (colored histograms). Golored circles mark the position of the bump 
for each population. Gontinuous lines superimposed on the histograms are the best-fitting 
Gaussians. 
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among the bumps is nearly constant in F438W, F606W, and F814W but increases by a 
factor of ~2 in F275W, as shown in Fig. [161 We also note that the order of the RGB-bump 
brightness is different in F336W, where population E still has extreme values of mFssew.bump, 
but the RGB bumps of both populations G and D are fainter than that of population B. 



Fig. 15.— Left panel: Reproduction of the mF 8 i 4 W vs. Ap 275 w,F 8 i 4 W diagram for RGB stars 
of Fig. 121 Middle panel: Zoom of left-panel diagram around the bump of population-E, where 
RGB-E stars are marked with blue dots. Right panel: Histogram of the F814W magnitude 
distribution. Grey error bars are Poisson errors while the red dashed-dotted line is the 
assumed luminosity function for the RGB (see text for details). 

To investigate the signihcance of the RGB bumps of populations B-E we have performed 
the analysis illustrated in Fig. [151 for population E. The left panel of Fig. [151 reproduces the 
verticalized mF 8 i 4 w vs. Ap 275 w F 814 W diagram of RGB stars shown in Fig. [21 while the middle 
panel is a zoom of the region around the bump of population E. The histogram in the right 
panel of Fig. [151 reproduces the F814W magnitude distribution for RGB-E stars (blue dots 
in the middle panel) already plotted in Fig. [HI The red line is the best £t straight line for 
the luminosity function of the upper-RGB and has been obtained by excluding stars with 
14.7 < mpsiAW < 15.1 to avoid the contamination from the bump. Grey error bars are 
Poisson errors calculated as the square root of the number of stars used to derive each point 
in the histogram minus one. 

Then we have simulated 1,000 verticalized GMDs for RGB stars by assuming same 
magnitude distribution as predicted by the red line, such that each synthetic diagram has 
the same number of RGB stars as in the observed GMD. We have applied to each GMD 
the same procedure to derive the F814W magnitude distribution as for real stars. For each 
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simulation we have calculated the difference between the area of the histogram derived from 
observations and the simulated one. 

/\ Avpn y^l4.8<mF8i4W<15.1 yyobserved simulated 

and find that AArea is signihcantly larger than zero in all the simulations. We extended the 
analysis to the bump of populations B, C, and D and find similar results thus concluding 
that the overdensities of stars resulting from the analysis illustrated in Fig. [Tl]are very likely 
the RGB bumps of populations B-E. 


In an attempt to interpret observations on the RGB bumps in the four populations we 
compare them with models. Specihcally, we have used two different sets of isochrones from 
BaSTI (Pietrinferni et ah 2004, 2006, 2009) and from Ventura et al. (1998, 2009, hereafter 
Roma models). 

We have estimated the luminosity difference between the RGB bump of helium-enhanced 
(Y=0.30, 0.35, 0.37) and helium-normal (Y=0.248) synthetic GMDs by using the same pro¬ 
cedure described above for real stars and compared these magnitude differences with ob¬ 
servations. Synthetic GMDs have been generated by using ASs and isochrones from both 
BaSTI and Roma models. We limited our study to the F814W filter which is only marginally 
affected by variations in light element abundances at hxed Z. We interpolated the value of 
Y which matches the observed and assumed the corresponding value of Y as the 

best estimate of the helium abundance of populations G, D, and E. 

The AY estimated from the difference in luminosity of the bump using both BaSTI and 
Roma models are listed in Table [5l We assumed that the stellar populations are coeval, as 
suggested by Piotto et al. (2007) who have revealed a narrow SGB in the analyzed mpsuw 
vs. mF 475 w — ''^F 8 i 4 W GMD, and used two different values for the absolute ages: 11.5 Gyr 
and 10.0 Gyr. We find that both population D and E are helium-enhanced with respect to 
population B, while population B and G share almost the same helium content, in analogy 
with what was inferred from multiple RGBs and MSs in Sects. 13.11 and 13.21 


6. Comparison among the AY from the different methods 

Figure fTTI shows the difference between the AY estimated from the RGB and the MS 
(upper panel), the bump location and the RGB (middle panel), and the bump location and 
the MS (lower panel). We note that there is a small discrepancy between the AY from the 
MS and the RGB for populations D and E, thought within <2a. The AY from the RGB 
is in fair agreement with the value calculated using the bump luminosity (within 1 a). The 
comparison of the middle panel of Fig. [T7| would favour a younger isochrones. 
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Fig. 16.— Yellow, cyan, and blue dots indicate the observed magnitude difference between 
the RGB-bump magnitude of populations C, D, and E (wix bump’ where X=F275W, F336W, 
F438W, F606W, and F814W), respectively, and the RGB-bump magnitude of population B 

(”^X,bump)- 


X 

^X.bumD 

^'X.bump 

™X.bump 

™X.bump 

F275W 

19.95±0.05 

19.94±0.02 

19.77±0.05 

19.59±0.05 

F336W 

17.75±0.03 

17.83±0.06 

17.78±0.03 

17.72±0.04 

F438W 

17.35±0.02 

17.36±0.02 

17.26±0.02 

17.18±0.04 

F606W 

15.99±0.02 

15.99±0.02 

15.91±0.02 

15.81±0.03 

F814W 

15.07±0.02 

15.09±0.02 

15.01±0.02 

14.91±0.03 

Model 

AVb 

AYc 

AVd 

AVb 

BaSTI, 10.0 Gyr 

0.000 

-0.010±0.008 

0.029±0.012 

0.103±0.025 

BaSTI, 11.5 Gyr 

0.000 

-0.008±0.008 

0.026±0.012 

0.098±0.025 

Roma, 10.0 Gyr 

0.000 

-0.009±0.013 

0.035±0.013 

0.100±0.023 

Roma, 11.5 Gyr 

0.000 

-0.007±0.013 

0.030±0.013 

0.085±0.022 


Table 5: Observed mx luminosity of RGB bump (X=F275W, F336W, F438W, F606W, 
F814W), and helium difference with respect to population B inferred from the F814W lumi¬ 
nosity of the bump. 
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Note that several authors have found that NGC 2808 has younger (10-15%) age with 
respect to the average age of intermediate-metallicity cluster (Rosenberg et ah 1999; De 
Angeli et ah 2005; VandenBerg et ah 2013; Milone et ah 2014). The comparison between 
the AY from the MS and the bump is somehow less satisfactory, though still within 1-2 a 
(depending on the adopted models). Again, a younger age seems to be favoured. A discussion 
on the origin of these discrepancies is beyond the purposes of the present paper and will be 
postponed to a future paper. 


7. The Horizontal Branch 

Among Galactic GGs, NGG2808 hosts one of the most extended HBs, which spans an 
extreme mpeoew — '^fsuw color range (L2=~0.9, Milone et ah 2014). The distribution of 
stars along the HB of this cluster exhibits three signihcant gaps, which separate four groups 
of HB stars. A red HB hosting about half of the total number of HB stars, and three distinct 
segments of blue HB stars (Sosin et ah 1997; Bedin et ah 2000; Piotto et ah 2007). The 
GMD shown in Fig. [T]con£rms this complex morphology. 

Spectroscopy of HB stars in NGG2808 further reveals that stars with different light- 
element abundance are distributed along different HB regions. Red-HB stars are Na-poor 
and 0-rich, while blue-HB stars are depleted in oxygen and enhanced in sodium (Gratton 
et al. 2011; Marino et ah 2014). 

The sodium distribution of red-HB stars is bimodal, with Na-rich stars having, on 
average, bluer B — V and U — V colors (Marino et al. 2014). The histogram of the [Na/Fe] 
distribution for red HB stars from Marino et al. (2014) is reproduced in the upper-left panel 
of Fig. [THl where we have colored red and blue the two stellar populations identified by these 
authors. For six HB stars in the HST held of view sodium abundances are available from 
Marino and collaborators. There is a clear anti-correlation between [Na/Fe] and mF 275 W 5 
with sodium-poor stars having also fainter luminosity in F275W as shown in the lower-left 
panel of Fig. [TSl 

To further investigate the connection between stellar populations with different sodium 
abundance and the red HB we combine optical and ultraviolet photometry. Since mF 275 w ~ 
mF 336 w and mF 336 w “ '^F 435 W colors are very efficient in separating stellar populations along 
the red HB of GGs, in the right panel of Fig. [13 we plot the mF 275 w vs. Gf 275 w,F 336 W,F 438 W 
Hess diagram for NGG2808. The fact that stars with different [Na/Fe] populate different 
regions of the mF 275 w vs. Gf 275 w,F 336 W,F 438 W diagram supports the conclusions by Gratton 
et al. (2011) and Marino et al. (2014) that the red HB of NGG2808 is not consistent with a 
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Fig. 17.— Upper panel: Difference between the AY estimated from the RGB and the MS 
for populations A, B, C, D, E. Middle panel: Difference between the AY coming from the 
bump location and the RGB. Lower panel: Difference between the AY coming from the 
bump location and the MS. We used two different isochrone sets and two different ages to 
estimate AY from the magnitude of the bump. 
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simple stellar population. 

The inset shows a zoom around the red HB. The distribution of stars along the HB 
is multimodal, with three main groups of stars clustered around Cf 275 W,F 336 W,F 438 W ~ 1-15, 
1.25, and 1.33 as highlighted by the histogram distribution of the pseudo-color Cf 275 W,F 336 W,F 438 W 
shown in the inset. An appropriate comparison with HB theoretical models is required to 
disentangle the effect of mass loss, evolved stars, and multiple stellar populations on the mor¬ 
phology of the red HB and to understand whether the three bumps correspond to distinct 
populations, but this is beyond the purposes of the present paper. 



-0.2 0 0.2 0.4 0 0.5 1 1.5 

[Na/Fe] F275W,F336W,F438W 

Fig. 18.— Upper-left panel: histogram of the distribution of [Na/Fe] for red HB 
stars (Marino et al. 2014). Lower-left panel: mF 275 w vs. [Na/Fe]. Right panel: mF 275 w 
vs. Cf 275 W,F 336 W,F 438 W Hess-diagram. The Hess diagram in the inset inset is a zoom around 
the red HB and the histogram of the distribution in Cf 275 W,F 336 W,F 438 W for red-HB stars is 
also shown. Na-rich and Na-poor stars dehned by Marino and collaborators are colored blue 
and red, respectively. 
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8. Multiple populations along the AGB 

Figure [19] illustrates our analysis of the AGB of NGC2808. The upper panels show a 
collection of mF 8 i 4 W vs. mx —'n7.F8i4W GMDs zoomed around the AGB (A=F275W, F336W, 
F438W, F606W). The 51 stars with ruFsuw > 13.5 that, on the basis of their position in 
these GMDs, are probable AGBs have been marked with colored symbols. 

The lower-left panel shows the mF 438 W vs. Gf 275 W,F 336 W,F 438 W diagram, where AGB stars 
are distributed along three distinct sequences. This feature is a signature of multiple stellar 
populations along the AGB. We dehned three groups of AGBi, AGBn, and AGBm stars 
that include 25, 11, and 15 stars, colored red, aqua, and magenta, respectively. These colors 
are used consistently in Fig. [191 AGBm is bluer than the remaining AGB stars in all the 
GMDs, while AGBu are slightly bluer than AGBi in all the GMDs apart from the mF 8 i 4 W 
vs. mF 336 w — fn-psiAW GMD where these two groups of AGB stars are almost overlapping. 
The mF 275 w vs. mF 275 w ~'^F 336 W GMD in the lower-right panel of Fig. [T^ shows that AGBm 
are brighter than the other AGB stars in F275W and that AGBu are, on average, brighter 
than AGBi. 

Theoretical models predict that hot HB stars would undergo a transition to an extended 
blueward nose excursion and exhibit bluer colors than the progeny of cold HB stars when 
reaching the AGB (e.g. Gingold 1976). It is tempting to speculate that the group of AGBm 
stars is the progeny of helium-rich HB, AGBi stars have primordial helium, and the AGBu 
belongs to a population with intermediate composition. The AGBi stars host 49±11% of the 
total number of AGB stars in agreement, within the large error bar, with the total fraction of 
the three helium-poorer RGB-A, RGB-B, RGB-G stars, which include half of the RGB stars 
of NGG 2808. The fraction of AGBm and AGBm stars are 22±6%, and 29±8%, respectively. 
These numbers only vaguely resemble the fraction of RGB-D (~31%) and RGB-E (~19%), 
though we admit there is some arbitrariness in selecting the AGB members of the three 
groups. 

We conclude that further analysis, possibly based on the synergy between spectroscopy 
and photometry, is needed to connect the triple AGB with the multiple RGBs of NGG 2808. 
Noticeably, stars with extremely thin H-rich envelopes miss the AGB phase, and move to¬ 
wards the white dwarf cooling sequence. It is not possible to hrmly establish from the present 
dataset if these ‘AGB Manque’ stars are present in NGG 2808 or not. 
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Fig. 19.— Upper panels: mpsuw vs. mx—irapsuw CMDs of AGB and RGB stars in NGG 2808 
(A=F275W, F336W, F438W, F606W). Lower panels: rnp/^^g^^ vs. Gp 275 'w,F 336 W,f 438 W (left) 
and mF 275 w vs. mF 275 w — ''^F 336 W GMD (right). Red, aqua, and magenta symbols represent 
the three groups of AGBi, AGBu, and AGBm stars dehned in the lower-left panel. 
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9. Discussion 

In this paper we have used multi-wavelength HST photometry to investigate multiple 
stellar populations in the GC NGC2808 as part of the Hubble Space Telescope UV Legacy 
Survey of Galactic GCs project (Piotto et ah 2015). 

Our basic results can be summarized as follows. 

• We have identihed hve distinct stellar groups along the RGB of NGG2808, namely 
A, B, G, D, and E, which contain 5.8±0.5%, 17.4±0.9%, 26.4±1.2%, 31.3±1.3%, and 
19.1±1.0% of the total number of RGB stars with 12.25 < mpsuw < 17.70, respec¬ 
tively. The hve stellar populations have been also detected along the MS althougth 
the separation between MS-A, MS-B, and MS-G is less evident than in the case of the 
RGB. We have found that the red MS discovered by Piotto et ah (2007) is composed of 
populations A, B, and G, while their middle and the blue MS correspond to population 
D and E, respectively. 

• We have exploited high-resolution spectroscopy from literature to infer the abundance 
of Na, O, Al, and Mg for the hve stellar populations. First of all, we have identihed 
the RGB-A-E stars for which chemical abundances from high-resolution spectroscopy 
are available. Specihcally, Garretta et al. (2006) and Garretta (2014) have found large 
star-to-star variations of [Na/Fe], [0/Fe], [Al,Fe], and [Mg/Fe] and identihed three 
groups of 0-normal, 0-poor, and 0-super-poor stars. We have matched the sample by 
Garretta and collaborators with our multi-wavelength photometry and found that 32 
of their stars belong to populations B, G, D, and E as dehned in this paper. 

Using the spectroscopic information, we found that population B has solar sodium-to- 
iron abundance ratio and is enhanced in oxygen ([O/Fe]~0.3 dex). Population G is 
enhanced in sodium ([Na/Fe]~0.2) dex and slightly depleted in oxygen by ~ —0.1 dex 
with respect to population B. Populations D and E are both sodium rich and oxygen 
poor and have [Na/Fe]~0.4 and [Na/Fe]~0.8, and [0/Fe]~ —0.4 and [0/Fe]~ —0.7, re¬ 
spectively. Unfortunately, no population-A RGB stars have spectroscopic information. 
Abundances of magnesium and aluminum are available only for three population-B 
and two population-G stars. All of them are distributed around [Mg/Fe]~0.4 and 
[Al/Fe]~0.1 with the two population-G stars being, on average, slightly enhanced in 
[Al/Fe] by ~0.1 dex. A larger sample is mandatory to establish if such a difference in 
Al is signihcant or not. 

• We have inferred the content of helium, G, N, O for the hve populations of NGG2808. 
To do this, we have followed the method by Milone et al. (2012b) and compared the 
observed colors with predictions from synthetic spectra. We found that populations 
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D and E are enhanced in helinm by ~0.06 and respectively, with respect to 

the populations B and C which share almost the same helium abundance. This he¬ 
lium difference follows from the assumption that population A and B have the same 
metallicity. From this assumption it also descend that population A has ~0.03 less 
helium than both population B and C. Planned spectroscopic observations in particu¬ 
lar of population A stars may allow us to test this assumption, and indicate whether 
the photometric differences between these two populations can be ascribed to iron and 
oxygen differences, rather than to helium. If we assume that population A has primor¬ 
dial helium (Y=0.246), our results indicate that population-E stars are highly helium 
enhanced up to V ~0.39. 

Note that it is the effect of helium-abundance variations on stellar temperatures that 
mostly causes optical and UV colors to change. Indeed stars with the same luminosity 
but different V have different effective temperature and gravity. On the contrary, 
helium has a marginal effect on the stellar atmosphere. (Sbordone et al. 2011). 

The comparison of observed and synthetic colors allow us to also estimate the average 
abundance of C, N, and O for each stellar population. We found that populations 
C, D, and E are enhanced in nitrogen by ~0.5, ~0.6, and ~0.8 dex with respect 
to population B, while population A has slight lower nitrogen (A[N/Fe]~0.1) than 
population B. Both oxygen and carbon anti-correlate with nitrogen. 

• We have detected the RGB bump of populations B, C, D, and E. In visual filters, 
which are marginally affected by light element-variations, the bump of population B 
and C have almost the same luminosity, while the RGB bump of population D and E 
are ~0.07 and ~0.17 mag brighter. The comparison of the observed bump luminosity 
with theoretical models suggests that populations D and E are more helium rich than 
both population B and G by ~0.03 and ~0.06 dex, respectively. 

• We confirm that the HB of NGG2808 is multi-modal with four main HB segments. In 
addition, we confirm that the red HB is inconsistent with a simple stellar population 
as suggested by the fact that the two groups of sodium-rich and sodium-poor stars 
identified by Marino et al. (2014) populate different regions along the red HB. 

• The AGB hosts three main sequences that are well distinguishable in the mF 438 w 
vs. Gf 275 W,F 336 W,F 438 W diagram. This finding indicates that the AGB of NGG2808 
hosts multiple stellar populations. 

In conclusion, the most astonishing property of this cluster is certainly its extremely 
complex stellar populations, as illustrated in Figs. [T]and|3l We have distinguished five dis¬ 
crete stellar populations, but a closer look at the Hess diagram shown in Fig. [3] suggests 
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that reality may be even more complex. Indeed, all five clumps appear to have some in¬ 
ternal structure, as if each of them could further split in two components, or having an 
internal spread in its photometric properties which may result from a small spread in chem¬ 
ical composition. This is particularly evident for populations B and C. We notice that our 
photometric data have an accuracy of ~ 0.01 mag, sufficient to resolve individual clumps, 
but distinguishing between an internal spread or a multiplicity would require larger samples 
of stars. 

The clear discreteness of the five (main) populations is a fact that every scenario for 
the formation of GCs must be able to account for. As emphasized in Paper I, this kind 
of discreteness is indeed an ubiquitous property among GCs and can not be ignored. It 
suggests that star formation occurred in a sequence of discrete events interleaved by periods of 
inactivity while the chemical composition of the interstellar medium was changing. However, 
we postpone to a future paper of this collaboration a dedicated discussion as to whether the 
various proposed scenarios can comply with this and the other observational constraints 
illustrated in Paper I. 

It also goes beyond the scope of this paper to try to identify in which temporal sequence 
the various multiple populations may have been generated, but we need to at least try to 
identify the first generation of this cluster. Population B is oxygen rich and sodium poor 
(see Fig. 8) and also is rather populous (~ 17% of the total in our sample), so it is the 
obvious candidate for being the first generation. However, if this is the case, we have a 
problem with population A, which is redder than population B in both mF 275 w ~ ''^fsuw 
and mF 336 W —''^F 438 W and for which no spectroscopic abundances are available from Garretta 
et al. (2006) and Garretta (2014). We have argued above that if Population B and A would 
have the same metallicity, than the helium abundance of population A should be lower that 
that of population B by ~ 0.03. Since no physical mechanism is known that could deplete on 
such scale the helium abundance below its Big Bang value, one would be forced to consider 
population A as the first generation. This is quite unpalatable, as population A represents 
only ~ 6% of the whole population of the cluster. Even if A-l-B together are regarded as the 
first generation, still they make only ~ 23% of the whole population sampled by WFG3 at 
the center of NGG2808. However, as in the case of other clusters, the first generation may 
be less centrally concentrated than subsequent generations and a more extensive mapping 
of this cluster is required to measure the overall fractions of the various populations. 

The alternative is to relax the assumption of these two populations having the same 
metallicity. Like in other clusters (e.g. M22, NGG 1851, M2, NGG5286), a small fraction of 
the core-collapse supernovae from the first generation may have contaminated the interstellar 
medium while such stellar population was still in the making. An increase of [Fe/H] by 
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~ 0.1 — 0.2 dex, associated to a parallel increase in oxygen as expected from core collapse 
supernovae, could then account for the photometric differences between population A and B. 
Our group is already engaged in high-resolution spectroscopic observations with GIRAFFE 
at the Very Large Telescope and spectra of stars in the hve RGBs of NGG2808 are going 
to be obtained soon. In particular, the abundances of iron, carbon, nitrogen, oxygen and 
sodium will be measured for population A, hence assessing whether this interpretation is 
viable. 

Moreover it should be noticed that, although some authors have suggested that multi¬ 
ple sequences are associated to distinct generations of stars, the possibility that GGs have 
experienced multiple or prolonged events of star formation is still strongly debated. We 
refer to papers by Bastian et al. 2013; Gabrera Ziri et ah 2014, 2015; Niederhofer et ah 2014; 
D’Ercole et al. 2008, 2010; D’Antona et al. 2005; Renzini et al. 2008; Decressin et al. 2007; 
Denissenkov et al. 2015 and references therein for critical discussion and for various scenarios 
and interpretations of multiple stellar populations. 
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